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The general solution of the equation "ty xx + - ^ x 4- Kx n e* = is displayed in terms of simple tabu- 
lated functions. The existence and uniqueness of solutions of a simple boundary value problem are 
determined as a function of the parameter K. 



For a class of problems in the steady flow of viscous 
fluids with a viscosity depending exponentially on 
temperature, an ordinary, nonlinear, second-order, 
differential equation arises [L2] 1 ; viz, 



V* 



* + -V x + Kx n e* = Q\x>Q. 
x 



(i) 



written in 
equations: 



the form of a system of two first order 



Kx n e*=f 2 {<p^,x) 



<Px 



--£-1 

X 



(3) 



The same equation may arise in calculating the tem- 
perature distribution in a dielectric in an alternating 
field [3]. In this equation, ^ is an unknown function of 
x, while K and n are arbitrary real constants. The 
equation is of considerable mathematical interest as an 
example of that exceedingly small class of strongly non- 
linear equations whose general solution can be dis- 
played in terms of simple tabulated functions. 

The reader may verify by direct differentiation that 
the following array is indeed a collection of solutions 
of eq (1): 



¥ 



= log{ 



Kx n+2 cosh 2 (m log x + a) 



^ = log 
^ = log 



-2m 2 



fc n+2 sinh 2 (mlog* + a)j ' 
-2m 2 



Kx n+2 sin 2 {m log x + a) J ' 

r ^ 1 

Kx;* +2 (log * + a) 2 J ' 



V=mlogx + a,K = 0. 



, K > 0. (2a) 
(2b') 
K < (2b") 

(2b"') 
(2c) 



In these equations, m and a are two arbitrary real 
constants which may be thought of as constants of 
integration. 

It is important to show that this array contains all 
the real solutions of eq (1). One can get at this ques- 
tion through some well known properties of ordinary 
differential equations. Observe that eq (1) can be 



It is immediately clear that the partial derivatives of 
f\ and/2 are continuous in the region x > 0, and hence, 
by a well-known theorem [4] eq (1) satisfies a Lip- 
schitz condition in some neighborhood of each point in 
the region. It follows immediately [5] that there is 
exactly one solution of eq (1) with a value and slope 
assigned at any given value of x ¥" 0. We can thus be 
sure that the array of eq (2) contains all the solutions 
if we can find values of the constants of integration 
which determine a solution of preassigned value and 
preassigned derivative for arbitrary, nonnegative x. 
We notice that eq (1) may be transformed in the follow- 
ing way. Let ^* = ^-f ^ and y = x/xo, then eq (1) 
becomes 

^* y + ±V* + Kxo*+*e-**y*e** = 0. (4) 

This equation is of the same form as eq (1) with K 
replaced by KxQ n+2 e~^ Q , while "^o and xo are unspecified 
constants. We see, therefore, that it is sufficient to 
show that among the solutions of eq (1) which vanish 
at x=l, there is a solution whose derivative at x=l 
takes on any preassigned value. 

Consider eq (2a) evaluated at x=l. By the vanish- 



ing of ^, we have that m- 



— cosh a, and thus that 



Wxh 



-(/i + 2) ±2 



sinh a. 



1 Figures in brackets indicate the literature references at the end of this paper. 



It is clear that we may choose a so as to make (^j-^i 
take on any real finite value. 
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When K < 0, any one of the three forms of the solu- 
tion, 2b', 2b", and 2b'" is possible. By the vanishing of 



"*¥ at x = 1 we have respectively, m = z 



sinh cz, 



m 



K a -»- /2 

— sin a and a=±W — — 

# = 1 can then be seen to be, 



K 



The slope of ^ at 



(^xXr-i = - (» + 2) ± 2 x /— cosh a, a # 



A 



(*x)«-i =-(fi + 2)±2 V— cos a, a^O 



Mr 



=-(fi + 2)±2V-7T 



2 



for the three forms 2b', 2b", and 2b"', respectively. 
It is again clear that we may choose a in these equa- 
tions so as to make {9 x )x=i take on any real finite 
value. We see further that all three forms are 
necessary, indeed, that we must choose form 2b', 
2b" or 2b'" according to whether the absolute value of 
0Pjp)jt=i + (n + 2) is greater than, less than, or equal to 

K 



For the case K = Q, the solution must be of the form 
(2c), but since eq (1) is linear when K — 0, the existence 
and uniqueness of solutions is well known and needs 
no further discussion. 

We have thus shown that the array of forms of eq (2) 
contains all real solutions of eq (1) in the region 
x>0. 

Of practical interest is the question of existence and 
uniqueness of solutions of boundary value problems. 
We shall consider a boundary value problem in which 
values of ^ are given at two values of x and a solution 
of eq (1) is sought which is bounded and smooth in the 
interval. Assuming that two solutions exist, call 
them ^i and ^2, we define the difference ^1— ^2 as a 
new function co. We then perform the following inte- 
gration by parts, 



L 



Xl d 

CO— (X(x) x )dx = JtCOCO., 



dx 



1 r *\ 
) J-''o 



xo) x 2 dx, 



(5) 



where x and X\ are the values of x at the boundaries of 
the interval. We replace the integral on the left of 
this equation by using the fact that both % and ^ 2 
satisfy eq (1) to get 



f* 1 

K x nJf 

Jx 



l e^2 W (l - e M )dx = xcocd.) 



T 



xo) x 2 dx. (6) 



Now it is very easy to see that the integrand on the left 
is always negative. Furthermore, if we restrict our- 
selves to solutions with bounded derivatives it is clear 



that the first term on the right vanishes (since co van- 
ishes on the boundaries) and, since the integral on the 
right is always positive, the right hand side of this 
equation is always negative. Hence, ifK is negative, 
there is at most one solution to the boundary value 
problem. If K is greater than zero this uniqueness 
proof breaks down, and indeed we shall see that mul- 
tiple solutions may be possible. 

We will next consider the case when K < 0, and 
show that there always exists a solution of the boun- 
dary value problem. We first use eq (4) to redefine 
K in terms of the given boundary data so that the 
vanishing of the scaled ^ at x = 1 corresponds to the 
boundary condition at the first boundary. When 
this condition is satisfied eq (2b) becomes 



sinh 2 a 



^ = log< 



^ = log < 



^ = log - 



^w+2 sm h2 ( a - 



— — log x sinh a 



sin^ a 



log x sin a) 



K 



~2~ log x± 1 



(7a) 



(7b) 



(7c) 



If the second boundary condition is ^ = ^2 at x = x>2, 
(x 2 > 1), we shall show that there is always one form of 
eq (7) which will satisfy the condition and which will 
define a smooth bounded function in the interval. For 

convenience, we define <p = exp -(^ 2 + (/i + 2)log%2). 

Then ip is a number determined by the boundary data 
which can be expressed in terms of the constant of 
integration a. From eq (7) we see that there are five 
possible forms of if. From (7a) we get two forms 2 



<Pi 



sinh a 



sinh I a + \~^~~ l°g ^2 sinh a J 



,0<a<<x, 



<£>2— " 



sinh a 



, < a < a, 



sinh I a — \~ir~ log X2 sinh a j 
from (7b) we get only one independent form 3 



<£>3— ■ 



. sin a 

sin I a -f \^~7T~ log x>i sin a j 



,0<a<7r, 



2 The symbol a here represents the positive root of the equation 



sinh a I- K 

— ^~ = ' V~2~ gX2, 



3 The form with the negative sign is equivalent to this form with a replaced by a + n. 
On the other hand <^, 3 (27r — a) = <p 3 {a), hence it is only necessary to consider values of a 
between zero and it. 
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and from (7c) we get two more forms 

1 



<£>4 : 



<P5 = 



K, 



—r- log x 2 + 



K ] 



— logX 2 - 



and again <p can take on any value between zero and 
infinity, so that there is always a solution to the 
boundary value problem. 

When K > 0, we have seen that the uniqueness proof 
fails. In fact we cannot expect a unique solution of the 
boundary value problem. Upon fitting the boundary 
condition ^P = at x= 1, eq (2a) becomes 



3sh 2 a 



x" +2 cosh 2 (a+ a/— — log .% cosh a) I 



00 < a< °o 



We confine our considerations, at first, to the situa- 



tion when 



*:, 



— — logX2< 1. A straightforward study 



of the range of values of the various forms of <p re- 
veals the following scheme: 



< & < 



1 



1 + 



log #2 



< <px < 



] 



-y lo g X2 



<(p 2 < 



in which, for a suitable choice of a and of the form of <p, 
(p can take on any value between zero and infinity. 
It is clear, therefore, that there is always a solution to 
the boundary value problem. 

f—~K 
When a/— — log x 2 ^ 1 some of the solutions rep- 
resented by eq (7) are no longer suitable because they 
no longer correspond to smooth bounded functions in 
the interval. In particular, forms (f 2 and <£>s correspond 
to inadmissible solutions, and in form <£> 3 , a must be 
restricted to the range 0< a< it — /3, where /3 is the 
smallest positive root of the equation 



Again, we represent the second boundary condition by 
ty = ty 2 at x = x 2 , {x 2 > 1). As before, we define 

(p = exp - (^2 + (n + 2) log x 2 ), and find 

cosh a 

T7~ (8) 



cosh(a -f- a/— log x 2 cosh a). 

A study of this function shows it to be a smooth posi- 
tive function of a, vanishing as a approaches plus or 
minus infinity and with a single maximum at a = a, 
where a is the root of the equation 



K, 



tanh {a + aJ— log x 2 cosh a) z 



tanh 



(9) 



A, 



f \/-^-log x 2 sinh a 



In this case, therefore, there are two solutions to the 
boundary value problem when <p < <£W f , one solution 
when <p = <pmax and no solutions when <p > <p ma x- This 
non-existence of solutions of the boundary-value prob- 
lem can be associated with an instability phenomenon 
in physical problems of viscous heating [6]. 
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sin f3 1/ -K 

We then find the following scheme for the possible 
values of the various forms of <p: 



< tpi < ■ 



1 



1 + 



■ < <£>3 < 



log* 2 
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